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The effects of quantum fluctuations in fields confined by background configurations may 
be simply and transparently computed using the Green's function approach pioneered 
by Schwinger. Not only can total energies and surface forces be computed in this way, 
but local energy densities, and in general, all components of the vacuum expectation 
on : value of the energy-momentum tensor may be calculated. For simple geometries this 

approach may be carried out exactly, which yields insight into what happens in less 
(N ■ tractable situations. In this talk I will concentrate on the example of a scalar field in a 

' circular cylindrical delta-function background. This situation is quite similar to that of a 

spherical delta-function background. The local energy density in these cases diverges as 
■ the surface of the background is approached, but these divergences are integrable. The 

^5 total energy is finite in strong coupling, but in weak coupling a divergence occurs in third 

I * ' order. This universal feature is shown to reflect a divergence in the energy associated 

I with the surface, the integrated local energy density within the shell itself, which surface 

energy should be removable by a process of renormalization. 

' Keywords: Casimir energy, divergences, renormalization 



1. Casimir Energies for Spheres and Cylinders 



- 1—1 

X 

The calculation of Casimir self-energies of material objects has become controver- 
sial, 1 although these concerns are nearly as old as the subject itself. 2-4 Although 
it appears possible to extract unique self-energies, they may be overwhelmed by 
terms which become divergent for ideal geometries. 5,6 Our attitude is that these 
terms may be uniquely removed by a process of renormalization, and that even the 
divergences revealed by heat-kernel methods 7 ' 8 may be unambiguously isolated. 

Table [1] summarizes the state of our knowledge concerning total Casimir self- 
energies for different simple configurations. The first row of the table refers to the 
Casimir energy of a perfectly conducting shell, either spherical or cylindrical, sub- 
ject to electromagnetic fluctuations in the exterior and interior regions. The second 
row refers to the same results for a scalar field subject to Dirichlet boundary con- 
ditions on the surface. The remaining four rows describe small perturbations: Row 
3 describes what happens for electromagnetic fluctuations when the interior of the 
sphere or cylinder is a dielectric having a permittivity e differing slightly from the 
vacuum value of unity; Row 4 indicates the same when the speed of light is the same 



1 



February 2, 2008 1:9 



WSPC - Proceedings Trim Size: 9.75in x 6.5in main 



2 



inside and outside the object, where £ = (V — e)/ (e' + e) in terms of the permittivity 
inside (e) and outside (e') the object; Row 5 shows the effect for a perfect conduc- 
tor of a small cllipticity Se (± refers to a prolate or oblate spheroid, respectively); 
and Row 6 refers to a f5-function potential (semitransparent shell) of strength A, 
which will be described in this paper. In these four cases, what is shown in the 
table is the coefficient of the second-order term in the relevant small quantity. One 
of the ongoing challenges facing quantum field theorists attempting to understand 
the quantum vacuum is to understand the pattern of signs and zeroes manifested 
in the table. 



Table 1. Casimir energy (E) for a sphere and Casimir en- 
ergy per unit length (£) for a cylinder, both of radius a. The 
signs indicate repulsion or attraction, respectively. 



Type 


^Sphere 11 


^Cylindcra 2 


References 


EM 


+0.04618 


-0.01356 


9, 10 


D 


+0.002817 


+0.0006148 


11, 12 


(e-l? 


+0.004767 = tMc- 

15367T 





13, 14 


e 


+0.04974 = ,f- 





15, 16 


&e 2 


±0.0009 





17, 18 


A 2 


+0.009947 = 





19, 20 



In this talk, we will illustrate the ideas for the interesting case of a circular 
cylindrically symmetric annular potential. Most of the calculations will refer to a 
(5-function potential. 



2. Green's Function 

We consider a massless scalar field in a r5-cylinder background, 



2a 



S(r — a) 



(1) 



a being the radius of the "semitransparent" cylinder. We recall that the massive 
case was earlier considered by Scandurra. 21 Note that with this definition, A is 
dimensionlcss. The time-Fourier transform of the Green's function, 



satisfies 



G(x,x') = J ^e--(*-*')<?(r,r'), 

g(r,r')=S(r-r'). 



-V 2 - lj z + -Sir - a) 
a 



Adopting cylindrical coordinates, we write 



C7(r,r')=||v fc ( z - Z '> jr ±.e im ^-^g m (r,r';k) 



(2) 



(3) 



(4) 



m— — oo 
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where the reduced Green's function satisfies 

,2 



1 d d 2 m X 
-— r— + k H — H — 5{r - a) 

r dr dr r z a 



l (r,r / ;fc) = -5(r-r'), (5) 
r 

where n 2 = k 2 — uj 2 . Let us immediately make a Euclidean rotation, 

u - iC, (6) 

where £ is real, so k is always real and positive. Apart from the 8 functions, Eq. ([5]) 
is the modified Bessel equation. 



2.1. Reduced Green's function 

Because of the Wronskian satisfied by the modified Bessel functions, 
K m {x)I' m {x) - K' m {x)I m {x) = -, 



(7) 



we have the general solution to the Green's function equation ([5|) as long as r =/= a 
to be 

g m (r,r';k) = I^Kr^K^Kr^ + A(r')I m (Kr) + B(r')K m (Kr), (8) 

where A and B are arbitrary functions of r'. Now we incorporate the effect of the 
S function at r = a in the Green's function equation. It implies that g m must be 
continuous at r = a, while it has a discontinuous derivative, 



a-r9 m (r, r';k) 
dr 



r— a+ 



= Xg m (a,r';k), 



(9) 



from which we rather immediately deduce the form of the Green's function inside 
and outside the cylinder: 



r,r' < a: g m (r, r'\ k) = / m (Kr < )ivT m (Kr > ) 



r,r > a 



1 + XI m (K,a)K m (na) 
g m (r,r';k) = 7 m ( r tr < )iT m (Kr>) 



{nr)I m {nr'), (10a) 
K m (Kr)K m ( K r'). (10b) 



1 + XI m (m)K m (Ka)' 

Notice that in the limit A — > oo we recover the Dirichlet cylinder result, that is, 
that g m vanishes at r = a. 



3. Pressure and Energy 

The easiest way to calculate the total energy is to compute the pressure on the 
cylindrical walls due to the quantum fluctuations in the field. This may be computed, 
at the one-loop level, from the vacuum expectation value of the stress tensor, 



(T" u ) = {d^d' v - ^g^d x d'^j ~G(x,a/) 



£(d»d v -g» v d 2 )~G(x,x). (11) 
i 
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Here we have included the conformal parameter £, which is equal to 1/6 for the 
conformal stress tensor. The conformal term does not contribute to the radial- 
radial component of the stress tensor, however, because then only transverse and 
time derivatives act on G(x, x), which depends only on r. The discontinuity of the 
expectation value of the radial-radial component of the stress tensor is the pressure 
on the cylindrical wall: 

P (^rr)in (^rr)out 

1 t r*r* xk1 



1671-3 m ~^oo ' -°° ^-oo 1 + XI m (Ka)K m (Ka) 



16tt 3 

m— — OO 



/*°° r°° ti d 
yldk dt- — ln[l + \I m (Ka)K m ( K a)], (12) 



where we've again used the Wronskian Q. Regarding ka and £a as the two Cartesian 
components of a two-dimensional vector, with magnitude x = na = y/ k 2 a 2 + C 2 a 2 , 
we get the stress on the cylinder per unit length to be 



1 f°° °° d 

S = 27raP=-- — j / dxx 2 > — In [1 + XI m (x)K m (x)] 
Aita 6 In dx 



(13) 



implying the Dirichlet limit as A — > oo. By integrating S = —J^£, we obtain the 
energy per unit length 

1 f°° °° d 

f = ~w/ dxx2 ^ln[l + A/ m (z)AT m (z)]. (14) 



m— — oo 



This formal expression will be regulated, and evaluated in weak and strong coupling, 
in the following. 

3.1. Energy 

Alternatively, we may compute the energy directly from the general formula 22 

To evaluate the energy in this case, we need the indefinite integrals 

dyyl 2 m {y) = \ [(x 2 + m 2 )I 2 m {x) - x 2 I 12 ] , (16a) 

1 1 

dyyK 2 n (y) = -- [(x 2 + m 2 )K 2 n (x) - x 2 K' 2 ] . (16b) 



When we insert the above construction (|T0|) of the Green's function, and perform 
the integrals as indicated over the regions interior and exterior to the cylinder, we 
obtain 

2 00 poo poo -i j 

£ = ~ J2 / d W dk( 2 -—\n{l + \I m (x)K m (x)}- (17) 

°^ ™_ „J-oo J —oo x ax 
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Again we regard the two integrals as over Cartesian coordinates, and replace the 
integral measure by 



dC dkC, 2 = ir / dnn 6 . (18) 

i J-oo JO 

The result for the energy (fT4|) immediately follows. 



4. Weak-coupling Evaluation 

Suppose we regard A as a small parameter, so let us expand the energy (I14j) in 
powers of A. The first term is 

£ I ^ -K m ( X )I m ( X ). (19) 
The addition theorem for the modified Bessel functions is 

oo 

K (kP) = e m( -^">K m (kp)I m (kp'), p>p', (20) 

m— — oo 

where P = \J p 1 + p' 2 — 2pp' cos(0 — 4>'). If this is extrapolated to the limit p' = p 
we conclude that the sum of the Bessel functions appearing in £W is K (0), that 
is, a constant, so there is no first-order contribution to the energy, £W = 0. 

4.1. Regulated numerical evaluation of £W 

Given that the above argument evidently formally omits divergent terms, it may 
be more satisfactory to offer a regulated numerical evaluation of £^> . We can very 
efficiently do so using the uniform asymptotic expansions (m — * oo): 



where x = rnz, t = + £ 2 , and = -jj- The polynomials in i appearing here 

are generated by 

1 r* i _ 5 S 2 

uo(t) = 1, u fc (t) = -t 2 {\ - i 2 K-i(*) + j ds — - — u fc _i(«). (22) 

Thus the asymptotic behavior of the products of Bessel functions appearing in 
Eq. ([in]) is obtained from 



fc=i 
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The first three polynomials occurring here are 



n(t) 
r 2 {t) 
rs{t) 



-(l-6t 2 + 5< 4 ), 



16 



(7 - 148i 2 + 554i 4 - 708t 6 + 295t 8 ), 



— (36 - 1666i 2 + 13775i 4 - 44272i 6 
16 



+ 67162i 8 - 48510t 10 + 13475£ 12 ). 
We regulate the sum and integral by inserting an exponential cutoff, S 



A 



Awa 2 



E' 

m=0 



dxx 2 —I m (x)K m (x)e xS , 
ax 



(24a) 
(24b) 

(24c) 
0+: 

(25) 



where the prime on the summation sign means that the m = term is counted with 
one-half weight. We break up this expression into five parts, 

A 



8na' 



■(I + II + III + IV + V). 



(26) 



The first term is the m = contribution, suitably subtracted to make it convergent 
(so the convergence factor may be omitted), 



I : 



dx x 



dx 



I (x)K (x) 



1 



2VT" 



-1. 



The second term is the above subtraction, 



II 



dx x' 



( d 



1 



\dx yT + x 2 



-xS 



26 



(27) 



(28) 



as may be verified by breaking the integral in two parts at A, 1 « A « 1/(5. The 
third term is the sum over the mth Bessel function with the two leading asymptotic 
approximants in Eq. (|23p subtracted: 



°° poo J 

m=2 y / dxx 2 4- 



I m {x)K m {x) - — 1 + —(1 - 6t 2 + 5t 4 ) 
Zm \ am z 



0. 



(29) 

Numerically, each term in the sum seems to be zero to machine accuracy. This is 
verified by computing the higher-order terms in that expansion, in terms of the 
polynomials in Eq. 



t 



t z 



I m {x)K m {x) - — 1 + — (1 - 6< 2 + 5t 4 ) 
Zm 



4m J 



r2(t)~-r 2 (t) 



8m 2 
t 

4m 7 



rz{t) ~\ri{t)r 2 {t) + l -r\(t) 



(30) 



which terms are easily seen to integrate to zero. The fourth term is the leading 
subtraction which appeared in the third term: 



OO pc 

IV = y m I 

m=l Jo 



dz z 



dz 



t e 



-mzS 



(31) 
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If we first carry out the sum on m we obtain 



dzr 3 



1 



1 



f 



1 1 



(32) 



,o (1 + ^ 2 ) 3/2 sinh 2 zS/2 S 2 25 6' 

as verified by breaking up the integral. The final term, due to the subleading sub 
traction, if unregulated, is the form of infinity times zero: 



8 ^ m 

rn—1 

Here the sum on m gives 



dzz 2 -^-(t 3 -6t° +5t Y )e 
az 



7\ „ — 7nz5 



oo 1 



In (1 



-zS 



)> 



and so we can write 



V = — — /" du (1 - u) Q W - 2 -«( U 3 / 2 - 6u 5 / 2 + 5u 7 / 2 ) 
16 eta J 



a=0 



Adding together these five terms we obtain 



0, 



(33) 



(34) 



(35) 



(36) 



87ra 2 f5 2 

that is, the 1/5 and constant terms cancel. The remaining divergence may be inter- 
preted as an irrelevant constant, since 5 — r/a, t being regarded as a point-splitting 
parameter. This thus agrees with the result stated at the beginning of this section. 



4.2. A 2 terra 

We can proceed the same way to evaluate the second-order contribution to Eq. ([i"4l 

\ 2 r°° r\ °° 

£(2) = it^ dxx dx E 

m— — oo 



(37) 



By squaring the sum rule (|20|) . and again taking the formal singular limit p' — > p, 
we evaluate the sum over Bessel functions appearing here as 

OO /.27T 7 

X] I 2 m {x)Kl{x)= -^ 2 (2xsm^/2). (38) 



m— — oo 



Then changing the order of integration, we can write the second-order energy as 

i2 /-27T 

dzzK$(z), (39) 



A 



dip 



647r 2 a 2 7 sin 2 <p/2 7 
where the Bessel-function integral has the value 1/2. However, the integral over ip 
is divergent. We interpret this integral by adopting an analytic regularization based 
on the integral (Res > —1) 

<•*" — 20Fr(i±^) 



d<P (sin |) 



r i 



(40) 
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Taking the right-side of this equation to define the <p integral for all s, we conclude 
that the tp integral, and hence the second-order energy £^ 2 \ is zero. 

The vanishing of the energy in order A and A 2 may be given a quite rigorous 
derivation in the zeta-function approach to Casimir energies — See Ref. 20. 



4.2.1. Alternative numerical evaluation 



Again we provide a numerical approach to bolster our argument. Subtracting and 
adding the leading asymptotic behaviors, we now write the second-order energy as 

[z = x/m) 



g{2) 



A 2 



87ra 2 



dx x 



I%{x)Kl(x) 



1 



4(1 + x 2 ) 



1 ™ r°° 2}- s r z i l °° 3 

- lim ^ 'mr s I dz + 2/ dzz— 

2--o±i Jo 1 + ^ 2 Jo *±lf=i 



r k (t) 



-,2k 



dx : 



I 2 m {x)K 2 m {x) - 



m=l J ° 

The successive terms are evaluated as 
A 2 



4m 2 



m—l k—1 



fc=i 



K , , A , 1, 9 C(2)^7C(4) 3106) 
- (7 + ln4)- i ln2 7 r- — + — - 



+0.000864 + 0.000006 



87ra- 



(0.000000), 



(41) 



(42) 



where in the last term in the energy (|41[) only the m—l and 2 terms are significant. 
Therefore, we have demonstrated numerically that the energy in order A 2 is zero to 
an accuracy of better than 10 -6 . 



4.2.2. Exponential regulator 

The astute listener will note that we used a standard, but possibly questionable, 
analytic regularization in defining the second term in energy above. Alternatively, 
as in Sec. l4.ll we could insert there an exponential regulator in each integral of e~ x6 , 
with 8 to be taken to zero at the end of the calculation. For m^Oi becomes mz, 
and then the sum on m becomes 

oo 

E ^ - ^T—y ( 43 ) 

m—l 

Then when we carry out the integral over z we obtain for that term 



7T 1, „ 
— In 2tt. 

85 4 



(44) 
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Thus we obtain the same finite part as above, but in addition an explicitly divergent 
term 



£ ( 2 ) = ^_ f 45 ) 

Again, if we think of the cutoff in terms of a vanishing proper time r, S = r/a, 
this divergent term is proportional to 1/a, so the divergence in the energy goes like 
L/a, if L is the (very large) length of the cylinder. This is of the form of the shape 
divergence encountered in Ref. 14. 



4.3. Divergence in 0(A 3 ) 

Although the first two orders in A identically vanish, a divergence in the energy (|14| 
does occur in 0(A 3 ). 

m=— oo u 

That such a divergence does occur generically in third order was proved in Ref. 20, 
using heat-kernel techniques. As we shall see, this divergence entirely arises from 
the surface energy. 



5. Strong Coupling 

The strong-coupling limit of the energy J14|) , that is, the Casimir energy of a Dirich- 
let cylinder, 



e° = - 1 



8?ra 2 ^ 

in— — oo ' 



f°° d 

J dxx 2 —\nI m {x)K m {x), (47) 



was worked out to high accuracy by Gosdzinski and Romeo, 12 



D _ 0.000614794033 
t - — . (48) 



It was later redone with less accuracy by Nesterenko and Pirozhenko. 23 For com- 
pleteness, let us sketch the evaluation here. Again subtracting and adding the lead- 
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ing asymptotics, we find for the energy per unit length 
1 



£ D = - 



8ira 2 \ 

OO 

m— 1 

oo 



- 2 



dx x 



dx x 

d 

dx 



ln(2xI (x)K (x))-l ' 



In (2xI m (x)K m (x)) - In 



1 + x 2 

xt\ ln(t) 



f°° d f d 

2V' / dxx 2 — \n2x + 2 V / dxx 2 — In 
^ Jo dx ^ Jr. dx 

°° f°° d 
> dxx 2 — 
±^ n ./o & 



mj 2 to 2 



1 1 



41+i 2 



1 OO 

2 ^ 

m—0 



dx- 



1 + x 2 



(49) 



In the first two terms we have subtracted the leading asymptotic behavior so the 
resulting integrals are convergent. Those terms are restored in the fourth, fifth, and 
sixth terms. The most divergent part of the Bessel functions are removed by the 
insertion of 2x in the corresponding integral, and its removal in the third term. 
(Elsewhere, such terms have been referred to as "contact terms.") The terms in- 
volving Bessel functions are evaluated numerically, where it is observed that the 
asymptotic value of the summand (for large m) in the second term is I/32to 2 . The 
fourth term is evaluated by writing it as 



21im V 



,2-s 



m— 1 



dz- 



1 + Z< 



- 2C(-2) = - 



03) 
2tt 2 ' 



(50) 



while the same argument, as anticipated, shows that the third "contact" term is 
zcroQ The sixth term is 

I" 

2 



— lim 

2 s^o 



C(s) 



- = Tln27r. 

s 4 



(51) 



The fifth term is elementary. The result then is 

£ D = -^r (0.010963 - 0.0227032 + + 0.0304485 + 0.21875 - 0.229735) 
0.0006146 



(52) 



which agrees with Eq. 



to the fourth significant figure. 



5.1. Exponential regulator 

As in the weak-coupling calculation, it may seem more satisfactory to insert an 
exponential regulator rather than use analytic regularization. Now it is the third, 
fourth, and sixth terms in the above expression that must be treated. The latter is 



a This argument is a bit suspect, since the analytic continuation that defines the integrals has no 
common region of existence. Thus the argument in the following subsection may be preferable. 
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just the negative of the term (f43|) encountered in weak coupling. We can combine 
the third and fourth terms to give 

1_ 2 f°° dzz 3 d 2 1 



6 2 



S 2 dz 2 



1 



(53) 



The latter integral may be evaluated by writing it as an integral along the entire z 
axis, and closing the contour in the upper half plane, thereby encircling the poles 
at iS and at 2inir, where n is a positive integer. The residue theorem then gives for 
that integral 

-P~2^> (54) 
so once again, comparing with Eq. (|50[) . we obtain the same finite part as in Eq. (|52[) . 

In this way of proceeding, then, in addition to the finite part found before in 
Eq. (f52")l . we obtain divergent terms 

1 1 1 



div 64a 2 5 



8ira 2 5 2 ' 4a 2 £ 3 ' ^ 
which, with the previous interpretation for <5, implies terms in the energy propor- 
tional to L/a (shape), L (length), and aL (area), respectively, and are therefore 
renormalizable. Had a logarithmic divergence occurred (as does occur in weak cou- 
pling in 0(A 3 )) such a renormalization would be impossible. However, see below! 



6. Local Energy Density 

We compute the energy density from the stress tensor (fTTj) . or 



( T oo) = I (d°d ' + V • V) G(x,x') 
2i 



-V 2 G(x,x) 
i 



16tt 3 



m— — oo 



— + d r d r> )g{r,r') 



2£-d r rd r g(r, r) 



(56) 



We omit the free part of the Green's function (jTU|) . since that corresponds to the 
vacuum energy in the absence of the cylinder. When we insert the remainder of the 
Green's function, we obtain the following expression for the energy density outside 
the cylindrical shell: 



u(r) 



16tt 3 



/OO rOO 
d( dk T 
-oo J — OO A 



+ XI m (Ka)K m (Ka) 



2uj 2 + K 2 + ^) K 2 m { K r) + K 2 K' 2 ( K r) 



-2S,-^r^K 2 m {Kr) 
r or Or 



r > a. 



(57) 
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The factor in square brackets can be easily seen to be, from the modified Bessel 
equation, 

2^Kl{nr) + l^i^K^r). (58) 

For the interior region, r < a, we have the corresponding expression for the energy 
density with I m <-> K m . 



6.1. Total and surface energy 

We first need to verify that we recover the expression for the energy found before. 
So let us integrate the above expression over the region exterior of the cylinder, 
and the corresponding interior expression over the inside region. The second term 
in Eq. ([58]) is a total derivative, while the first may be integrated according to the 
integrals given in Eq. p^|) . In fact that term is exactly that evaluated above. The 
result is 

/(dr)u(r) = / dxx 2 — + XI m {x)K m {x)] 



The first term is the total energy (TH| . but what do we make of the second term? 
In strong coupling, it would represent a constant that should have no physical 
significance (a contact term — it is independent of a if we revert to the physical 
variable k as the integration variable) . 

In general, however, there is another contribution to the total energy, residing 
precisely on the singular surface. This surface energy is given in general by 22 ' 24-28 



l-4g 
2i 



<j> dS ■ VG(x,x') 



(60) 



which turns out to be the negative of the second term in J(dr) u(r) given in Eq. (|59[) . 
This is an example of the general theorem 

(fir) u(r) + £ = E, (61) 

that is, the total energy E is the sum of the integrated local energy density and the 
surface energy. A consequence of this theorem is that the total energy, unlike the 
local energy density, is independent of the conformal parameter £. 



6.2. Surface divergences 

We now turn to an examination of the behavior of the local energy density as r 
approaches a from outside the cylinder. To do this we use the uniform asymptotic 
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expansion (|2ip . Let us begin by considering the strong-coupling limit, a Dirich- 
let cylinder. If we stop with only the leading asymptotic behavior, we obtain the 
expression (z = nr/m) 



where 



cLk k 



oo 

E < 

m=— oo 



9 Tit , _ 9 7T 1 

-k 2 + 2(1-4£)k 2 

2to 2r7?i z 2 



(A 



X = -2 



77(2) - 77 (z^) 



(62) 



(63) 



and we have carried out the "angular" integral as in Eq. (|18p . Here we ignore the 
difference between r and a except in the exponent, and we now replace k by mz/a. 
Close to the surface, 



2 r - a 



X ~ T 



f r 

and we carry out the sum over m according to 

d 3 1 12 



3 t 4 r 



4^4 



c?X 3 X X 4 4 (r - a) 4 ' 



Then the energy density behaves, as r 

u(r) 



1 



16tt 2 (r - a) 



;(l-6£). 



(64) 



(65) 



(66) 



This is the universal surface divergence first discovered by Deutsch and Candelas. 2 
It therefore occurs, with precisely the same numerical coefficient, near a Dirich- 
let plate 19 or a Dirichlet sphere. 29 It is utterly without physical significance (in 
the absence of gravity), and may be eliminated with the conformal choice for the 
parameter £, £ = 1/6. 



6.3. Conformal surface divergence 

We will henceforth make this conformal choice. Then the leading divergence depends 
upon the curvature. This was also worked out by Deutsch and Candelas; 2 for the 
case of a cylinder, that result is 

u(r) ~ 9 —. — rrr , T — > a+, (67) 

exactly 1/2 that for a Dirichlet sphere of radius a. To get this result, we keep the 
1/m corrections in the uniform asymptotic expansion, and the next term in x' 

2r — a ( r — a \ 2 
X~7— + M— • (68) 
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6.4. Weak coupling 

Let us now expand the energy density (|57|) for small coupling, 



u(r) 



\ />00 />OC ^ ^ 

-— / dc dk ^hE(-^)"chch 

107T J — oo J —oo ™_ ^ „_n 



-K 2 + (1 - 4£) K 2 + 



^( K r) + (1-40« 2 ^M . (69) 



If we again use the leading uniform asymptotic expansions for the Bessel functions 
we obtain the expression for the leading behavior of the term of order A n , 



u (n) (r) 



— k—, r / ffizz > 



The sum on m is asymptotic to 



m 3 ~ n e 



-mx 



(3-n)! 



2(r - a) 



r — > a+, 



so the most singular behavior of the order A™ term is, as r — ► a+, 



tt (»)( r ) ~ (-A) f 



(3- w)!(l -6£) 
967r 2 r"(r - a) 4 "™' 



(71) 



(72) 



This is exactly the result found for the weak-coupling limit for a <5-sphere 29 and for 
a (5-plane, 22 so this is a universal result, without physical significance. It may be 
made to vanish by choosing the conformal value £ = 1/6. 



6.5. Conformal weak coupling 

With this conformal choice, once again we must expand to higher order. Besides 
the corrections noted in Sec. 16.31 we also need 



t = t(za/r) ~ t + (t- i 3 )- — -, r->a, (73) 



Then a quite simple calculation gives 



(«) 



(n-l)(n+2)r(3-n) 
( ~ A) 28807r 2 r»+i(r-a)3-' T ~* a+ ' (?4) 



which is analytically continued from the region 1 < Ren < 3. Remarkably, this 
is exactly one-half the result found in the same weak-coupling expansion for the 
leading conformal divergence outside a sphere. 29 Therefore, like the strong-coupling 
result, this limit is universal, depending on the sum of the principal curvatures of the 
interface. Note this vanishes for n = 1, so in every case this divergence is integrable. 
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7. Cylindrical Shell of Finite Thickness 

We now regard the shell (annulus) to have a finite thickness 6. We consider the 
potential 

A„t = ~<?<r{r), (75) 

where 

{0, r<a_, 
h,a-<r< a+, (76) 
0, a + < r. 

Here a± = a±5/2, and we set hS = 1. In the limit as (5 — > we recover the (^-function 
potential. As for the sphere 29 it is straightforward to find the Green's function for 
this potential. In fact, the result may be obtained from the reduced Green's function 
given in Ref. 29 by an evident substitution. Here, we content ourselves by stating 
the result for the Green's function in the region of the annulus, a_ < r, r' < a + : 

g m {r,r') = 7 m (/tV<)ii' m ( r tV>) + AI m (K'r)I m (n'r') 

+ B{I m (Kr)K m (Kr') + K m (n 'r)I m {n' V')] + CK m (K,'r)K m (K,'r'), 

(77) 

where k' = \J n 2 + Xh/a. The coefficients appearing here are 

A = -^\nl' m {Ka-)K m {K' a-) - k' ' I m {Ka-)K' m {n' a_)] 

x[KK' m (na + )K m (n' a+) - n ' K m {na + )K' m (K 'a+)\, (78a) 

B = ^[Kl' m (Ka-)I m (K'a-) - K'I m (Ka-)I' m (K'a-)} 

x [K,K' m {na + )K m (K a+) - k ' K m {na + )K' m (K ' a+)\, (78b) 

C = -^[Kl' m {Ka-)I m (K' a-) - k! I m {Ka-)I' m (n' a-)] 

x [KK' m (Ka + )I m (K'a+) - k' K m {Ka+)I' m {K' a+)\, (78c) 



where the denominator is 



5 = [Kl' m (Ka-)K m (K' a-) - k ' I m (na-)K' m (K 

X [ K ^m( Ka +) / m( K ' a +) - K '^m(Ka+)/^(K'a + )] 

- [K7^(/ta_)/ m (/t'a_) - K7 m (Ka_)/'„(K'a_)] 

x [K J fT r '„(Ka + ) J fT m (rt'a + ) - «;'ii' m (Ka + ) J fT r '„(K'a + )]. (79) 
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7.1. Energy within the shell 

The general expression for the energy density within the shell is given in terms of 
these coefficients by 



u(r) 



1 



*** jo 



dn k 



-k 2 + (1-40 



1 d d 



r dr dr 

WlM'r) + CK 2 m {K'r) + 2BK m (K'r)I m (K'r)}. (80) 



7.2. Leading surface divergence 

The above expressions are somewhat formidable. Therefore, to isolate the most di- 
vergent structure, we replace the Bessel functions by the leading uniform asymptotic 
behavior (f2Tj) . A simple calculation implies 



.4 
B 
C 



t+ + 1' + 

t + - 1' + i_ - 1'_ 
t+ + 1' t- + 1'_ 



,2m(ri'_-v'+) 



t_ - t'_ 
t- + t'_ 



(81a) 
(81b) 
(81c) 



where t + = t(z+), rj_ = rj(z'_), z'_ = n'ci-fm, etc. If we now insert this approxima- 
tion into the form for the energy density, we find 



u = (T uu ) 



8tt 



| 00 PC 

a + m=l Jo 



dz-\- z+t r 



t+ - 1 



t+ + 



t'J + 



,2,2 



■(1-80 + 



t- - t'_ 

tl + 1'_ 
4 



,2m(-r)' T +v'-) 



(1-4C) 



m z 



+ — + 1 + — +' 

2^2 b + «■+<'- <-- „2m(T}L-»7^) 



(82) 



If we are interested in the surface divergence as r approaches the outer radius 
a + from within the annulus, the dominant term comes from the first exponential 
factor only. Because we are considering the limit Xha «m ! ,w have 

A/ ' " " ' ' (83) 



t', w t+ I 1 
+ + 1 2m 2 



and we have 



Xh/a 



32ttVI ^ 

~r m— 1 



m/ d»«t(l -8^ + t 2 )e 2m «-0. 



(84) 
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The sum over m is carried out according to Eq. (fTTj) , or 



E 



me 



2m(ri' r -ri' + ) 



rtL 



2{r-a+) 



and the remaining integrals over z are elementary. The result is 

Aft 1 - 6£ 



96ir 2 a (r — &+) 2 



r -> o+, 
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(85) 



(86) 



the expected universal divergence of a scalar held near a surface of discontinuity, 30 
without significance, which may be eliminated by setting £ = 1/6. 



7.3. Surface energy 

Now we want to establish that the surface energy £ (|6"0j) is the same as the integrated 
local energy density in the annulus when the limit S — > is taken. To examine this 
limit, we consider Aft/a 3> k 2 ■ So we apply the uniform asymptotic expansion for the 
Bessel functions of k' only. We must keep the first two terms in powers of k <C k': 



, 2 I rn (Ka-)K m (na + ) , 
a ~ — k : — ; — , — smhm(rj_ — rf + ) 



mz'_z' +y Jt'_t' + 
1 ftT 1 

— \ 7r I ' m ( Ka ~) K rr t (na+) - — J -±I m (Ka-)K' m (Ka + ) 

Z j y -2> V 

x coshm(?7'_ — rf,). 



Because we are now regarding the shell as very thin, 

S 1 



where 



V- - V+ 



at' 



'Aha 



using the Wronskian ((7|) we get the denominator 

5 ~ — \[l + \I m (m)K m (na)}. 



Then we immediately find the interior coefficients 
.4 



B 

C 



7r \ / Afta Im(Ka)K m (Ka) c _ 2m ,/ 
2 1 + XI m (Ka)K m (m) 

2 1 + AI m (na)K m (na) ' 

1 



— VXha 
2tt 



I m {na)K m (Ka) ^ 2mn > 
1 + XI m (Ka)K m (na) 



(87) 
(88) 
(89) 

(90) 

(91a) 
(91b) 
(91c) 
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7.4. Identity of shell energy and surface energy 

We now insert this in the expression for the energy density (|80[) and keep only the 
largest terms, thereby neglecting k 2 relative to Xh/a. This gives a leading term 
proportional to h, which when multiplied by the area of the annulus 27ra<5 gives for 
the energy in the shell 



which is exactly the form of the surface energy 2; given by the negative of the second 
term in the integrated energy density ([59")) . 

7.5. Renormalizability of surface energy 

In particular, note that the term in l£ of order A 3 is, for the conformal value £ = 1/6, 
exactly equal to that term in the total energy £ in Eq. (|4^|) : 



This means that the divergence encountered in the global energy is exactly ac- 
counted for by the divergence in the surface energy, which would seem to provide 
strong evidence in favor of the renormalizablity of that divergence. 

8. Conclusion 

The work reported here and in Refs. 20,29 represents a significant advance in under- 
standing the divergence structure of Casimir self-energies. We have shown that the 
surface energy of a (5-function shell potential is in fact the integrated local energy 
density contained within the shell when the shell is given a finite thickness. That 
surface energy contains the entire third-order divergence in the total Casimir energy. 
The local Casimir energy diverges as the shell is approached, but that divergence is 
integrable, so it yields a finite contribution to the total energy. The identification of 
the divergent part of the total energy with that associated with the surface strongly 
suggests that this divergence can be absorbed in a renormalization of parameters 
describing the background potential. 

Challenges yet remain. This renormalization procedure needs to be made precise. 
Further, we must make more progress in understanding the sign (and for cylindrical 
geometries, the vanishing) of the total Casimir self-energy. And, of course, we must 
understand the implications of surface divergences on the coupling to gravity. Work 
is proceeding in all these directions. 
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(92) 



g(3) = £(3), 



(93) 
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